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ABSTRACT 

Methods for solving Maxwell’s equations are integral part of optical metrology and computational lithography 
setups. Applications require accurate geometrical resolution, high numerical accuracy and/or low computation 
times. We present a finite-element based electromagnetic field solver relying on unstructured 3D meshes and 
adaptive hp-refinement. We apply the method for simulating light scattering off arrays of high aspect-ratio 
nano-posts and FinFETs. 

Keywords: Scatterometry, optical metrology, computational metrology, computational lithography, 3D rigorous 
electromagnetic field simulations, finite-element methods, hp-FEM 

1. INTRODUCTION 

Optical metrology can be used to detect features sizes on a sub-nanometer level. In the semiconductor industry 
it is used in process control and in mask quality control for pushing performance of DUV and EUV lithography. El 
Numerical modelling is an important part of optical metrology setups in this field: measurement results are 
compared to simulation results of a parameterized model in order to quantitatively determine dimensions of 
the measured sample. With increasing complexity and decreasing feature sizes, the need for accurate optical 

Corresponding author: S. Burger 
URL: http://www.jcmwave.com 
URL: http://www.zib.de 



Figure 1. FinFET: Left: Images of the geometry layout in two cross-sections. Center: Images of parts of the tetrahedral 
mesh in different viewing directions. Right: Computed reflection spectra for S- and P-polarized light. 
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Figure 2. Left: Schematics of the geometry of the investigated scatterometric target (unit cell of a 2D-periodic grating of 
nano-posts with circular cross section, not to scale). Different parameters of the model are indicated (critical dimension, 
CD, at different heights, hi, h 2 , pitches px, py). Center / right: Reflection spectra of S/P-polarized light, -Rs/p,pi for 
numerical discretization parameters p = 


metrology methods for complex 3D shapes is increasingPI^ This triggers also a need for efficient numerical 
methods for computational metrology. 

The challenge for electromagnetic field (EMF) solvers (Maxwell solvers) is typically efficiency (i.e., to achieve 
highly accurate results at low computation times). Finite-element methods (FFM) allow for high efficiency due 
to accurate geometrical modelling, adaptive meshing strategies, and higher-order convergence. In simulation 
tasks requiring high accuracy FFM can outperform other rigorous simulation methods.^^^ 

We develop and investigate finite-element methods for electromagnetic field simulations. In previous contri¬ 
butions these have been applied to various setups in optical metrologyIn this context, also 3D structures 
have been investigated.^^^lti^ Here we discuss methods for further performance improvements, especially for 
efficient simulation of 3D devices with complex geometries. This is reached by using hp-finite elements on un¬ 
structured, tetrahedral and prismatoidal meshes. Figure [1] shows the model, mesh, and simulation result of a 
typical investigated sample. 

This paper is structured as follows: The background of the numerical method is presented in Section [5] 
The method is validated by presenting simulation results for two different examples related to current critical 
dimension (CD) metrology requirements Section[3]presents simulations of nano-post arrays with high aspect- 
ratio. Section m presents simulations of arrays of fin field-effect transistors (FinFET). 

2. HP FINITE ELEMENT METHOD 

In the following the background of the finite element method is summarized.^ Light scattering off nanoscopic 
structures on scatterometry samples is modeled by the linear Maxwell’s equations in frequency domain.l^l^ 
From these a single equation for the electric field E can be derived: 

curl p,“^curl E — w^eE = iwJ. (I) 

where e and p, are the permittivity and permeability tensor, oj is the time-harmonic frequency of the electromag¬ 
netic field, and the electric current J is source of an electromagnetic field. The domain of interest is separated 
into an infinite exterior Dext which hosts the given incident field and the scattered field, and an interior Hint where 
the total field is computed. Electromagnetic waves incident from the exterior to the interior at the boundaries 
between both domains are added to the right hand side of Eq. O- For numerical simulations the infinite exterior 
is treated using transparent boundary conditions (using the perfectly matched layer method, PML). 

For a FFM discretization, Eq. © is first transformed into a weak formulation, i.e., it is tested with a vectorial 
function (j) and integrated over which yields: 
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Figure 3. Spectral dependence of numerical discretization error ARs/p_p for S/P polarization, obtained for various 
horizontal FEM degrees p, with fixed Pz = 5. 
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For compact notation, the forms a{4>, E) and f{ip) are introduced, and the function space H (curl) is defined!^ 
The weak form of Maxwell’s equations then reads: 

Find E G H (curl) such that: 

a(,^,E) = /(<(.), V,^GH(curl). (3) 


A finite element discretization of Maxwell’s equations restricts the formulation (|31) to a finite-dimensional 
subspace 14 with dim 14 = N < oo: 

Find E/i G 14 such that: 

a{(l>h,'Eh) = fiPh). ^(ph&Vh. (4) 


Next, a basis {ipi, ..., ^pN} of 14 is constructed, and the electric field is expanded using the basis elements: 

N 

^h = Y: Ciipi. The variational problem (Eq. [T]) is then tested with all elements of the basis which gives a linear 
system of equations: 

N 

'^a{<pj,<pi)ei = f{ipj), Vj = l,...,Af. (5) 

i=l 

The matrix Aji = a{ipj,(pi) is sparse and can be decomposed with efficient sparse LU solvers to obtain the 
unknown expansion coefficients of the electric field. 

The basis {ipi, ■ ■ ■, is constructed using elements pi (also called ansatz functions) which are polynomial 
functions of order p, and which are defined on a single patch of the spatial discretization of the geometry (mesh) 
only. For the results presented in this paper, we attribute elements of different polynomial order p to different 
patches of the mesh.^ In regions where the mesh is very fine due to required geometry resolution (very thin 
layers or other fine details of the geometry) a lower polynomial order p can be chosen than in regions where 
the mesh is coarser. The method to distribute different orders p to the different patches relies on estimating 
errors on the different patches, making use of informations on geometry, meshing, material properties and source 
fields. This yields a basis {pi,..., pjv} which is well adapted to the problem and does not spend too much 
computational effort in regions where it is not required. 

We demonstrate this so called hp-FEM method for two different applications: an array of very elongated 
posts and a transistor geometry with nanometer features (FinFET). For the array of nano-posts, due to the 
elongated geometry, best performance is reached when the mesh consists of elongated elements (in this case 


















prismatoidal elements) and when different polynomial degrees p are used in the dimension of elongation and 
in the orthogonal dimensions. For the transistor, an unstructured tetrahedral mesh discretizes the geometry, 
where the typical dimensions of the mesh elements, h, can vary over about two orders of magnitude. In this case 
different polynomial degrees p are used for different tetrahedral elements. 

3. SIMULATION OF LIGHT SCATTERING OFF SILICON NANO-POSTS WITH 

HIGH ASPECT-RATIO 

The model investigated in this section corresponds to a periodic array of Silicon nano-post on a silicon sub¬ 
strate. The geometrical setup is described schematically in Fig. [5] (left), and the parameter configuration for the 
simulations is defined in Tabled] 


Px = Py 

32 nm 

CDi / CD2 / CDa 

16 nm / 22 nm / 24 nm 

hi jhi 

300 nm / 400 nm 

d! cf 

30 deg / Odeg 


Table 1. Parameter settings for the Si nano-post array (compare Fig. [2|). 


In the model the structure is illuminated from above (superspace with refractive index n = 1) with S- and 
P-polarized plane waves, at an angle 0 to the surface normal, and rotation angle (j). For simulating a wavelength 
spectrum, independent computations of the time-harmonic model for 100 wavelengths are performed. The 
optical material parameters of Silicon at each wavelength are obtained from tabulated data.^^ In our setup, up 
to 80 independent computations are processed in parallel on the computation threads of a standard workstation. 
Typical computation times for a single computation at low to medium accuracy levels range from below 1 sec to 
few seconds. 

Fig. m (center, right) shows spectra of reflected S- and P-polarized light, Rs and Rp. The displayed spectra 
are obtained for the same physical setting, however, at different numerical resolutions p, i.e., Rs,p. In this case, 
we always use the same spatial mesh where the geometry is discretized using prismatoidal elements which are 
elongated in z-direction (surface normal). The numerical discretization parameter p here corresponds to the 
polynomial order of the finite-element ansatz functions pi (cf., Sec.[2|) in the x — y-plane. The polynomial order 
of the ansatz functions in z-direction in this case is chosen as Pz = 5. For the absolute values on a linear 
scale, as displayed in Fig. [S] differences in the computed spectra for p > 2 can hardly be detected. Therefore, 
in Figure [3| the differences of the computed spectra to the spectra computed at highest numerical resolution. 



Figure 4. Left: Convergence of the average error ARp with horizontal finite element degree p, for S- and P-polarized 
incident light. Right: Convergence of the average error ARp with (geometrical) discretization parameter Nc (number of 
segments of the polygon defining the cross-section of the nano-post). 



















^R-S,p = |R-s.p ~ Rs,p= 5 |j are displayed on a logarithmic scale. AR is also termed numerical discretization error. 
As can be expected the numerical discretization error decreases exponentially with p in the whole spectral range. 
We define the average discretization error as ARp = ARp/N, where the summation is performed over all 
N spectral points. Please note that several different conventions are used for defining numerical discretization 
errors in this context. The average error obtained from the data in Fig. [3] is displayed in Fig. |4] (left). In the 
investigated parameter regime, exponential convergence to very high accuracy levels is observed. 

All results displayed in Figures [5] to S] (left) are obtained on the same discretization of the geometry (essen¬ 
tially on the same prismatoidal mesh). In order to verify accuracy of the discrete geometry we have performed a 
convergence study for a meshing parameter: For meshing the nano-post geometry the mesh generator automat¬ 
ically discretizes the circular cross-section of the nano-post with a polygon, with a given number of segments, 
Nc- We have computed reflectivity spectra similar as in Fig. [5] for different segment numbers Nc- Here we have 
defined the average discretization error as AR^r^ with AR^r^ = jR^r^ — RAr^= 36 |. As can be seen 

from Fig. 0] (right) the results converge very well with geometry discretization parameter Nc- Results with an 
accuracy AR^r^ < 0.2% are obtained for W > 10. 

4. SIMULATION OF LIGHT SCATTERING OFF FINFETS 

The model investigated in this section corresponds to a periodic array of FinFET structures of the 22 nm (and 
smaller) technology nodes. Bunday et aP point out that with the launch of such small structures on integrated 
circuits complex 3D architectures have become a cruicial driver for down-scaling, and that this also implies 
additional needs for metrology. 

The investigated geometrical setup is schematically shown in Fig. [T] (left). All dimensions of the device follow 
Fig. 15 and Table 10 of Bunday et aP (22 nm node: Fin pitch: 44 nm. Gate pitch: 88 nm, Fin width: 12.7 nm. 
Fin SWA: 89.5 deg. Fin height: 40 nm. Fin undercut: 2.1 nm, Gate width: 40 nm. Gate SWA: 89.8 deg. Gate 
height: 95 nm, Gate undercut: 2 nm, SiN thickness: 5 nm, High-k layer thickness: 2 nm, TiN thickness: 7 nm, 
BOX thickness: 200nm). Additional parameters are roundings of the fin (rp: 3nm) and of the gate (re: 8nm) 
top edges. The optical material parameters in the investigated spectral range are again obtained from tabulated 
dateP (Si, SiN, SiOz, TiN, TiOz). 

Fig. [1] (center) shows parts of tetrahedral meshes discretizing the FinFET geometry. Meshing is performed 
with an automatic mesh-generator,which is a part of the finite-element package. The unstructured mesh 
allows to accurately resolve fine geometrical features like sub-nm undercuts, sidewall-angles and corner-roundings 
with only few additional mesh elements, facilitating accurate geometry modelling. When high refractive index 
contrasts between the different involved materials are present, such accurate geometry resolution is essential for 




Figure 5. Hp-adaptive FinFET simulations: Left: Fraction of mesh elements on which ansatz functions <pi with polyomial 
order p are used, as a function of the numerical accuracy parameter pprec. Right: Convergence of average numerical error 
ARpprg^ with accuracy parameter pprec- 




















Figure 6. Tetrahedral meshes of the FinFET geometry with different settings of corner roundings: Left: No corner 
rounding, rp = 0, re = 0. Center: rp = Inm, re = 8nm. Right: rp = 4nm, re = 8nm. 


precise approximation of the electro-magnetic field distributions. This enables accurate light scattering results 
in the far field, as inaccurate near-held resolution is generally not ’’smoothed out” in the far held. Computed 
rehection spectra for S- and P-polarized incident light at oblique angle of incidence are displayed in Fig. [T] (right). 

The shortest dimensions of the edges of the tetrahedra in the meshes displayed in Fig. [1] (center) are com¬ 
parable to the smallest geometrical features, i.e., typically smaller than Inm. However, for best performance, 
in regions of larger geometrical structures (e.g., in the gate, or in the buried oxide layer, etc.) mesh element 
dimensions are rather scaled with some fractions of the wavelength of light, i.e., several orders of magnitude 
larger than the dimensions of the smallest mesh elements. 



Figure 7. Left: Reflection spectra for P-polarized light, without corner rounding (rp = 0, re = 0) and with corner 
rounding (rp = 3nm, re = 8nm). Right: Difference between the two reflection spectra, dRp = Rp(rp = 0,rG = 
0) — i?p(rp = 3nm, re = 8nm). 
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Figure 8. Spectral sensitivity of reflection of S (left) and P (right) polarized light for changes of fin corner rounding radius 
rp around the value of rp = 3nm. 


When asssembling the finite-element matrix Aji = a{(pj,(pi) (cf. Eq. [^, a-priori error-estimation is used to 
choose the individual polynomial order p of the set of ansatz functions (pi on each individual mesh element. A 
parameter pprec controls error-estimation. Figure [S] (left) shows how with different settings of pprec the fraction 
of mesh elements where different polynomial orders p are used is changing. E.g., for a setting of pprec = 10“^ 
on about 60% of the mesh elements, ansatz functions pi with polyomial order p = 1 are used, on about 35% 
of the mesh elements, ansatz functions with p = 2 are used, and on the remaining elements, ansatz functions 
with p > 3 are used. With decreasing pprec these percentages shift towards higher p for the individual patches. 
Therefore, the dimension of Aji increases with decreasing pprec, numerical discretization errors are expected to 
decrease with pprec, and computational costs increase with pprec, however, only moderately.l^ Figure [5] (right) 
shows convergence of the computed spectrum with pprec^ very high accuracy levels are obtained. 

Unstructured meshes and finite elements with adaptive (hp) choice of polynomial order p allows to obtain 
fast and accurate results. This allows to analyze and reconstruct fine geometry features for complex, multi-scale 
geometries. To demonstrate this we investigate the influence of corner rounding on scattering spectra. Sunday 
et aP note fin corner rounding as a parameter which does influence FinFET performance. Figure 0 compares 
spectra for FinFET geometries with and without fin and gate corner roundings. Figure [S] shows the changes of 
S and P reflection spectra with fin corner rounding radius rp. In these simulations, gate rounding was fixed to 
tg = 8 nm, and fin rounding was varied by 1 nm around a rounding radius of rp = 3 nm. The investigations can 
be used to design the scatterometric measurement setup for best sensitivity for specific parameters of interest .1^ 
We note that sensitivities can also be computed directly at high accuracies.^^ Values of dR in Fig. |S] are of the 
order of up to few times 10“^ for changes of the rounding radius of 1 nm. This suggests that comparable or 
higher numerical accuarcy ranges should be used in reconstructions of such parameters. 

5. CONCLUSION 

A finite-element method using hp-adaptivity on 3D prismatoidal and tetrahedral meshes has been demonstrated. 
Convergence to highly accurate results over a wide spectral range has been observed for examples related to CD 
metrology challenges. The method allows to efficiently compute scatterometric signals and sensitivities. 
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